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We analyze the amplification due to so-called superradiance from the scattering of pulses off 
■ rotating black holes as a numerical time evolution problem. We consider the "worst possible case" 

' of scalar field pulses for which superradiance effects yield amplifications < 1%. We show that this 

— , , small effect can be isolated by numerically evolving quasi-monochromatic, modulated pulses with 

_ ■ a recently developed Teukolsky code. The results show that it is possible to study superradiance 

in the time domain, but only if the initial data is carefully tuned. This illustrates the intrinsic 
difficulties of detecting superradiance in more general evolution scenarios. 

m ; 

I. INTRODUCTION 

7— I . 

' In the last few years, we have been involved in the development of a numerical code for the time evolution of 
0^ , perturbations of rotating black holes based on the Teukolsky equation |l]J^]. There are several motivations for this 
work. One is the desire to revisit problems that have previously (mainly in the 1970s) only been approached in the 
frequency domain. That is, our goal is to explore the effects of the rotation of the black hole from a "time-evolution" 
point of view. More importantly, our ultimate goal is to provide a framework that will, once we understand how to 
construct astrophysically relevant initial data H, be used to extend the close-limit approximation of head-on black 
Q\ . hole collisions to the case of inspiral black hole mergers. The working premise in head-on, close- limit collisions Q| 
""^ ' is that the merger can be viewed as perturbations of non-rotating black holes. In contrast, an inspiral close-limit 
Q-( approximation requires perturbations about a rotating black hole. 

Our Teukolsky code project took us first to study the dynamics of scalar fields in the Kerr geometry jij. This work 
mainly concerned the late-time, power-law behaviour of a scalar perturbation. The second installment concerned 
gravitational perturbations and discussed the full dynamical response of a black hole to an external perturbation, 
namely the quasinormal mode ringing and the subsequent late-time tails. In Ref. we also dealt briefly with super- 
radiance: The anticipated amplification as certain wavelengths are scattered by the black hole. However, although the 
results we obtained indicated the presence of superradiance , we feel that our previous analysis was not completely 
satisfactory. Hence, the goal of this short paper is to return to the issue of superradiance in a setting that yields 
unequivocal evidence for the superradiance phenomenon. 

The direct approach to measure superradiance from the time evolution of perturbations of rotating black holes 
is to compute the energy flux going "down the hole" . For perturbative fields that posses well-defined stress-energy 
tensors (e.g. scalar and electromagnetic fields), it is possible to construct such a conserved energy flux ||. The case of 
gravitational perturbations is not that simple ]5||| . For this reason, we will concentrate our analysis on the "simple" 
case of scalar perturbations. The price to pay is that superradiant effects in this case are < 1% Q], thus requiring a 
highly accurate evolution code. 

For scalar perturbations, the Teukolsky equation in Boyer-Lindquist coordinates reads 
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Above, M is the mass of the black hole, a is its angular momentum per unit mass and A = r 2 — 2Mr + a 2 . The two 
horizons of the black hole follow from A = 0, and correspond to r± = M ± v M 2 — a 2 . Reference to the azimuthal 
angle tp has been removed by assuming that the perturbation has a harmonic dependence e" nip . 
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II. SUPERRADIANCE IN THE FREQUENCY DOMAIN 



In the standard approach to solve the Teukolsky equation, one proceeds via separation of variables. For our present 
purposes, it is sufficient to note that this essentially corresponds to assuming that i) the time-dependence of the 
perturbation is accounted for via Fourier transformation, and ii) there exists a suitable set of angular function that 
can be used to separate the coordinates r and 9. In the case of scalar perturbations, the angular functions turn out 
to be standard spheroidal wave- functions J8|. Knowing this, we assume a representation (for each given integer m) 
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where it should be noted that the angular functions depend explicitly on the frequency uj; that is, they are intrinsically 
time-dependent functions. After a separation of variables of the form given by (|^), the problem reduces to a single 
ordinary differential equation for Ri m (r,u>). This equation can be written as 
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where K = (r 2 + a 2 )ui — am, G — rA/ (r 2 



2 ) 2 , and the tortoise coordinate r» is defined from <ir* = (r 2 +a 2 )/Adr. 
The variable E is the angular separation constant. In the limiting case a — > 0, it reduces to 1(1 + 1), and for nonzero 
a it can be obtained from a power series in aui ||. It should be noted that E is real valued for real frequencies. 
The physical solution to (||) is defined by the asymptotic behaviour 
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where uj + = a/2Mr + is the angular velocity of the event horizon. T and 1Z denote the transmission and reflection 
coefficients, respectively. Given the prescribed asymptotic behaviour, together with that for the complex conjugate of 
Rim and the fact that two linearly independent solutions to (|J) must lead to a constant Wronskian, it is not difficult 
to show that 



(1 — muj + /uj)T =1 — 7?. 
From this result, it is evident that superradiance (1Z > 1) occurs if 



uj < muj+ 



2Mr+ 
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Alternatively, one can deduce that energy can be extracted from the black hole immediately from the boundary 
condition (|J). If u> < muj+, the solution oc e - l (^- rnuJ +) r * ) which behaves as "ingoing" into the horizon according to 
a local observer, will in fact correspond to waves coming out of the hole according to an observer at infinity. That is, 
for superradiant frequencies, one would expect to find energy flowing out from the horizon, cf. JTcj ], Superradiance is 
the wave analogue of the standard Penrose process, and its existence implies that it would in principle be possible to 
mine a rotating black hole for some of its rotational energy. This may seem exciting, but it is very unlikely that this 
effect will play a relevant role in any reasonable astrophysical scenario. Nonetheless, it is an interesting effect that 
deserves a close theoretical investigation. 

In Fig. [I], we show a sample of results for the reflection coefficient in the case when I = m = 2. These results were 
obtained by a straightforward integration of (Q) and subsequent extraction of 1Z. The maximum amplification in this 
case is close to 0.2 %. Similar results were obtained by Teukolsky and Press more than 25 years ago They also 

considered electromagnetic waves and gravitational perturbations and found that the maximum amplification is 0.3% 
for scalar waves, 4.4% for electromagnetic waves and as large as 138% for gravitational waves. 

Given the results in Fig. [l], it is worth pointing out that they agree with the standard conclusions regarding the 



apparent "size" of a rotating black hole as seen by different observers. It is well-known (see for example 11 ) that 



the black hole will appear larger to a particle moving around it in a retrograde orbit than to a particle in a prograde 
orbit. This is illustrated by the fact that the unstable circular photon orbit (at r = 3M in the non-rotating case) 
is located at r = AM for a retrograde photon, while it lies at r = M for a prograde photon. The results in Fig. [j] 
illustrate the same effect: In our case, we have prograde motion when oj/m is positive and retrograde motion when 
ui/m is negative. The data in Fig. |l| correspond to m = 2, and the enhanced reflection for positive frequencies as 
a — > M has the effect that the black hole "looks smaller" to such waves. Conversely, the slightly decreased reflection 
for negative frequencies leads to the black hole appearing "larger" as a — > M. 
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III. SUPERRADIANCE IN THE TIME DOMAIN 



From the discussion in the previous section, it should be clear what must be done to construct initial data sets 
that would yield superradiance with our Teukolsky code. Following an idea introduced in Ref. 0, we construct 
superradiant initial data by setting up a pulse containing mainly frequencies in the interval < ui < muj + (from 
now on we will assume that m is positive). Obviously a non-superradiant pulse would be one whose main frequency 
content is outside this window. Furthermore, it makes the analysis easier and better suited to comparisons with 
frequency domain calculations if the pulse is "almost monochromatic." To achieve this, we use as initial data an 
ingoing Gaussian pulse modulated by a monochromatic wave. Assuming that this initial pulse is centered far away 
from the black hole at = r Q and that the modulation frequency is er, we set at t = 

<j> of e -(r,-r a +t) 2 /b 2 -ia(r,-r +t) _ (j^ 

Thus, the corresponding power spectrum is P{u>) = P max e _ ( a,-<7 ) b / 4 . The modulated pulse leads to a Gaussian 
frequency distribution that peaks at a frequency u) = a. Given this initial data, and the fact that there is no 
frequency dispersion in the Teukolsky equation, we ought to be able to detect superradiance in our evolutions if 
< a < mu>-\-. Conversely, if a > muj + or a < we should not find any amplification in the scattered waves. Finally, 
it is not enough for the peak of the power spectrum of the pulse to lie within the superradiant frequency window. To 
maximize the effect, we need also to minimize the "frequency overlap" of the initial pulse into the non-superradiant 
regime. To accomplish this, we have at our disposal the parameter b that governs the width of the pulse. If for 
instance we want P(m uj+) / P max — e, we should use 



b= h/^m. (8) 

muj + — a 

Before showing that the pulse (0) can indeed be used to probe superradiance in the time domain, a few comments 
on our Teukolsky code are needed. The code was described in detail in Ref. jjj, but there is one specific issue that is 
important for the present study that has not yet been discussed. To avoid numerical difficulties (cf. the description 
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in Of]), we replace the azimuthal angle ip with the "ingoing Kerr-coordinate" , which is defined by 

<P = <P+ I ^dr . (9) 

The transformation between the solution $ to the standard Teukolsky equation and the solution $ from our numerical 
code in terms of <p is 

(10) 
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This is more or less obvious, but it is important to notice a few things. First of all, the replacement ip — ► <p 
changes the symmetry of the equations. While the original Teukolsky equation (l|) is symmetric under the change 
(m, ip) — > (—m, —p), the equation we evolve with our code for <p is not. That is, while evolutions for the same Gaussian 
pulse (unmodulated!) should lead to the same emerging scalar waves for ±m in the original case, this will not happen 
when we use <p. To ensure that the anticipated symmetries are present in our results, we have constructed initial 
data in Boyer-Lindquist coordinates and then used the transformation ( |Io|) to get data in the tp coordinate system. 
This is especially important in the case of modulated Gaussians since the idea that we could enhance superradiance 
in the scattered wave by centering the Gaussian at a non-zero frequency was based on an analysis in Boyer-Lindquist 
coordinates, and the result of the experiment will depend on a careful tuning of the initial pulse. 

Armed with the above conclusions, we are prepared to discuss our numerical results. Because we are interested in 
unveiling the amplification due to superradiance, we will focus on the energy flux through various surfaces surrounding 
the black hole. 

Given a spacetime with a time Killing vector t a (like the Kerr geometry) and a perturbation with a well-defined 
stress-energy tensor T a b, it is possible to define || a conserved energy flux vector T a bt b . The flux of energy across a 
3-dimensional time-like hypersurface with unit normal r a is then given by 

dE = T ab t a r b dS , (11) 

where dS is the 3-surface element of the hypersurface. For a massless scalar field, 
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with over-bars denoting complex conjugation. For simplicity, we monitor the flux of energy through r = constant 
surfaces in Boyer-Lindquist coordinates. This assumption together with r a r a — 1 yield r a = ±(0, \fK/ p, 0, 0), where 
(as before) A = r 2 — 2Mr + a 2 and p 2 = r 2 + a 2 cos 2 9. Furthermore, the time Killing vector in this case reads 
t a = (1, 0, 0, 0), and the surface element is given explicitly by 



dS = \/-g {3) sin9d6dipdt = y/Ap $m9d9dipdt . (13) 
By collecting the above results and noticing that t a r a — 0, it is not difficult to show that 

dE = ±-(d r $dt$ + dr&dt$)ABmed9d(pdt . (14) 

Finally, integration over ip yields 

dE = ±tt (<9 r $a t $ + a r $9 t #)A smOdOdt , (15) 
which can be rewritten in terms of the tortoise coordinate as 

dE = ±n{d rt $d t <f> + 9 r ,$5 t $)(r 2 + a 2 ) sin (9 d9 dt . (16) 

In our simulations, the energy flux is monitored through two surfaces located at r* = ±20M. The outer surface is 
well away from the black hole while the inner one is reasonably close to the event horizon. The scattering of waves 
by the curved spacetime should be strongest, i.e. most of the interesting physics should have its origin, in the region 
included between these surfaces. 

In Figs. H and ||, we show the results of the superradiance "experiment" . The displayed data are for two qualitatively 
different situations. Both datasets correspond to evolutions with m = 2 and a black hole rotation parameter a = 
0.99M. In both cases the pulse was initially centered around r a — 125M, and the angular distribution of the initial 
data was chosen to be the standard spherical harmonic Y 2 (9, (f). The first case (see Fig. ||) shows a situation where 
one would expect to see superradiance. We have chosen the modulation frequency of the impinging Gaussian (see the 
previous section) such that a — 0.8 mid+, and the width of the Gaussian corresponds to e = 0.01. That is, the pulse 
has its main support (in the frequency domain) in the superradiant regime. The second case (Fig. |3|) corresponds to 
a Gaussian with the same width but now centered around a = -0.8mw+. 

As is obvious from Fig. l], the scattering of the two pulses we have chosen should be quite different. And not 
surprisingly, we find that this is indeed the case. In the superradiant case shown in Fig. 0, all the initial energy is 
reflected by the black hole. Superradiance is distinguished in two ways. By monitoring the energy flowing across 
the surface at r* = 20M, we see that the reflected energy is slightly amplified after scattering (cf. the upper panel 
of Fig. ||). In this specific case, the amplification corresponds to 0.14%. It should be compared to the maximum 
single frequency amplification of 0.187% for a = 0.99M, deduced from the data in Fig. [j] (and also the 0.3% found by 
Teukolsky and Press ||). That we are seeing superradiance is also clear from the fact that energy flows out through 
the surface at r* = — 20 M (cf. the lower panel of Fig. ||). The total energy flowing out through the inner surface 
corresponds to a superradiant amplification of 0.11%, in reasonable agreement with the result deduced at the outer 
surface. 

The non-superradiant results, obtained by modulating the Gaussian with a frequency of opposite sign to that used 
in the superradiant case, are in clear contrast to the superradiant ones. In Fig. |^, there is certainly no amplification 
of the reflected wave. In fact, as can be seen from the upper panel of Fig. @, the infalling pulse is almost entirely 
swallowed by the black hole. That there would be very little reflection in this case could, of course, be anticipated by 
comparing our chosen Gaussian pulse to the data in Fig. [l| 



IV. CONCLUDING REMARKS 



We have designed a numerical experiment that clearly exhibits the presence of superradiance phenomenon when 
waves of a certain character are scattered by a rotating black hole. Superradiant effects have previously only been 
studied in the frequency domain, and one conclusion that can be drawn from the present work is that superradiance 
is perhaps best approached in that way. True, we have managed to extract the amplification due to superradiance 
in the "worst possible case" of scalar waves, when the superradiant amplification is expected to be considerably less 
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than 1%, but this was mainly due to having at our disposal a conserved flux. More than anything else this is direct 
evidence of the precision of our evolution code to solve the Teukolsky equation (|J|] . 

The initial data we used to isolate the tiny effect due to superradiance is required to be "almost monochromatic" and 
thus artificial. Our numerical experiment shows that superradiance can play a role in evolutions when the scattered 
pulse has support only in a restricted frequency domain. This is undoubtedly an interesting illustration, but what 
about superradiance in more general cases? It seems to us that the effect is easiest to isolate if one monitors different 
frequencies separately, i.e. works in the frequency domain. The main reason for this is that an amplification of a 
reflected signal with increasing a is not in itself an indication of superradiance. The results shown in Fig. [I] indicate 
that one would generally expect enhanced reflection of prograde moving waves as a — > M. This effect is likely to 
overwhelm the actual "amplification" of certain superradiant frequencies in an evolution of general initial data. This 
is certainly true for scalar waves, and since the maximum amplification of impinging electromagnetic waves is only a 
few percent [^) the conclusion should hold also in that case. However, the possibility that superradiance may play a 
distinctive role in an "astrophysical" evolution for gravitational waves cannot be ruled out. For gravitational waves, 
the amplitude of certain frequencies should be amplified by more than a factor of two || . A detailed study of that case 
could provide interesting results, but the present results for scalar field prompts us to proceed with caution. We have 
learned that the initial data requires careful tuning in order that superradiance be observed. For general initial data, 
absorption of the non-superradiant frequencies will typically make the amplification due to superradiance difficult to 
distinguish. Moreover, superradiance should not be confused with the competing effect that the "size" of the black 
hole changes with the rate of rotation. As we have seen, this effect will generally lead to a much enhanced reflection 
of prograde waves which may confuse an attempt to distinguish superradiance. 
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FIG. 1. Reflection coefficient (7V) as a function of the frequency of the wave (cu M) for different values of the angular momen- 
tum parameter (o) with I = m — 2. Superradiance is present in the interval < w M < ma/2r + , with r + — M + V ' M 2 — a 2 . 
The bottom panel is a close-up of this superradiance regime. The maximum observed superradiant amplification is ~ 0.187%. 
As a — * M, there is a clearly enhanced reflection of prograde waves (ui > 0) while the overall reflection of retrograde waves 
(w < 0) decreases somewhat. 
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FIG. 2. An example of a superradiant evolution. We show results corresponding to a modulated Gaussian pulse (with 
a = 0.8mcj + , m = 2, a — 0.99M and e = 0.01, see the main text for more details). We monitor the integrated energy flux 
through two surfaces, one located at r* = 20M and the second at r» = — 20M. (The normals of these surfaces are chosen such 
that energy flowing inwards into-the-hole across the outer surface is positive.) Superradiance manifests itself in two ways. At 
the outer surface (the upper panel), we see a total amplification of ~ 0.14% in the reflected energy. As is clear from the inset 
in the upper panel, the reflected energy is larger than what initially fell onto the black hole. At the inner surface (the lower 
panel), we find that energy mainly flows out of the horizon. The integrated flux of this energy corresponds to ~ 0.11% of the 
total energy falling onto the hole, i.e. agrees well with the amount that amplifies the reflected waves at the outer surface. 
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FIG. 3. An example of a non-superradiant evolution. The data are similar to that described in Fig. 2, but here the 
modulation frequency of the initial Gaussian is a — -0.8mu + , and there is no sign of superradiance. 



8 



